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Abstract 

There are now many nodal superconductors in heavy fermion (HF) systems, 
charge conjugated organic metals, high Tc cuprates and ruthenates. On the 
other hand only few of them have well established A(k). We present here a 
study of the angular dependent thermal conductivity in the vortex state of 
some of the nodal superconductors. We hope it will help to identify the nodal 
directions in A(k) of UPd2Al3 , UNi2Al3, UBeia and URu2Si2. 
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Unconventional superconductivity moved to center stage in recent years |T[. After a long 
controversy d-wave symmetry of both hole- doped and electron- doped high cuprates has 
now been established Among HF superconductors the symmetry of A(k) has been 

identified only for UPts |QJ^ and very recently CoCelus 0. As to the organic supercon- 
ductors the symmetry of k-(ET)2Cu(NCS)2 has also been established recently through the 
angular dependent thermal conductivity in the vortex state Indeed we have shown 

that the angular dependent thermal conductivity provides a unique window to access the 
nodal directions in A(k) |]8|-|T2|. For example from study of angular dependent thermal 
conductivity in Sr2Ru04 in a magnetic field within the a-b plane, we conclude that the 



nodal lines do not lie within the a-b plane but along the c- axis possibly at ckz=±^ p!3| , P 
Recently a number of authors proposed that the nodal lines in UPd2Al3 are not in the a-b 
plane (ck^=0) but somewhere along the c axis |T5|-p!8| (ck^ 7^0). The presence of the nodal 
lines is clearly seen from the low temperature thermal conductivity data [|l5i|l9| . Also since 
UPd2Al3 has the hexagonal symmetry, the nodal lines should lie along the c axis. The 
most obvious possibilities and their associated angular variation of the thermal conductiv- 
ity will be the focus of this work. It is generally assumed that UPd2Al3 is a spin singlet 



superconductor since it shows a Knight shift reduction below T^ p0|-p2|. Its magnitude is 
however difficult to interpret due to contributions of the AF ordered local moments. The 
excitations of the latter have been identified to mediate the electron pairing in UPd2Al3 



18| . Theoretical analysis of this mechanism shows that due to the crystalline electric 



field anisotropy a nondegenerate odd parity state should be most favorable. It also has a 
reduced spin susceptibility below Tc and an associated Knight shift. Therefore we include 
one odd parity nodal gap function with ck2=0 in our discussion. We want to stress that the 
theory we will develop here will also help to determine the nodal lines in UBei3, URu2Si2 
and UNi2Al3 which has recently been proposed to exhibit odd parity spin triplet pairing 
21. 



In the following we compute the quasiparticle density of states and the thermal con- 
ductivity in a vortex state at low temperature (T-C A) in a magnetic field within the a-c 
plane. Also we limit ourselves to quasi two dimensional systems where the Fermi surface is 
approximated by a corrugated cylinder. When the ratio of Fermi velocities f3=Vc/^a is not 
too small, the angular dependent thermal conductivity will tell the nodal position along the 
k^- axis. 



To calculate the quasiparticle density of states we focus ourselves to nodal superconduc- 
tors where the nodes lie at some values of ckz=±Xo Also for the usefulness of our results 
tt=(vc/va)^ should not be too small (say a >0.1). A magnetic field H is applied in the a-c 
plane with an angle 6 from the c axis. The supercurrent around vortices is circulating in 



the plane perpendicular to H. Following Volovik |25], we can handle the effect of the super- 



current within semiclassical approximation. Then the residual density of states is given by 



9 = Re{ ^ \ (1) 

^ /(Co-ix)2 + f2/ 
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where x= |v ■ q|/A and f= A(k)/A, Co = \imi^^Qlm{uj/A), u is the renormahzed 
frequency and k the quasiparticle wave vector. Also v • q is the Doppler shift with v the 
quasiparticle velocity and 2q the pair momentum, the brackets in the equation above denote 
the average over the vortex lattice and the Fermi surface. As to f we limit ourselves to three 
simple cases with odd and even parity gap functions respectively which are appropriate 
for a Fermi surface model with cylindrical symmetry 1) f=sm{ckz), 2) f=cos(cfcz) and 3) 
f=cos(2cfcz). Then for these three cases the average is worked out and we obtain 

2 / 2 X 

g = -(Coln(^^^=) +a;tan-^(— ; 



~ (x) + Co((ln(^)) - 1) for {x) » Co (2) 
^Coln(A) + for (x)«Co 

We call (x) 3>Co and (x) ^Cq the superclean limit and the clean limit respectively 
||9|JTl[|. Here (...) means the spatial average over the vortex lattice. Now let us assume that 
the quasiparticle relaxation is due to impurity scattering in the unitarity limit. Then Co is 
determined by 

C, ^ f (3) 

Which is solved as 

Co = (f |l"(2y'H|))i(l - :^{x')) for (i) « Co (4) 
Finally the residual density of states is given by 

= (7o(l + -^(x2)) for (x) <Co (5) 



and 



yu 2A ° ^2A^ ^ 

(x) = -VaV^ne) (6) 

TT 

1 9A 

{x') = -vl{eH)F^{e)H- 



with 
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FIG. 1. Angular functions Fa{0) and Fc{9) in the clean limit and isotropic case (a=l) for 
different node line positions xo=ckz. 

I{9) = {cos^ 9 + asin^ 6)-^- f dcf) 

n Jo 

[cos^ 9 + sin^ ^(sin^ (p + a sin^ xo) 

sinxo cos0sm(26')]2 (7) 

Fc{9) = (cos^ 9 + a sin^ 6^)^ [cos^ 9 + {^ + a sin^ xo) sin^ 

where a=(vc/va)^ and Xo= 0, | and | for f=sin(cfc^), cos{ckz) and cos{2ckz) respectively. 
The 9 dependence of Fc{9) and Gc{9)=P{9) is shown as for x=0, f and ^ in the isotropic 
case a=l in Figs. 1 and 2 (dashed lines). Boundary values for Xo=0 are given by Fc(|) = 
^y/a = |(^) and Gc(f ) = Similar curves for the anisotropic case with a=0.6 are 

shown in Figs. 4 and 5 respectively. As noted elsewhere the 9 dependent residual density 
of states is accessible to the specific heat, the superfiuid density and the spin susceptibility 
PJTT[]. For example 

Cs/inT = g 

Xs/xn = g (8) 
pJps{H = 0) = l-g 

where C^, Xs and ps are the electronic specific heat, spin susceptibility and superfiuid 
density in the vortex state respectively. Furthermore 7iv is the Sommerfeld coefficient and 
Xn is the Pauli susceptibility in the normal state. 

Explicitly the field dependent specific heat is given by 

Cs/inT = T — 1{9) 

vr A 
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FIG. 2. Angular functions Ga{0)= l{9)i{9) and Gc{0)= 1(6')^ in the superclean limit and 
isotropic case (a=l) for different node line positions xo=ck2. 

for the superclean limit (9) 
for the clean limit 

In the latter the Cs{9) angle dependence is given directly by Fc(6') shown in Figs. 1 and 
4. For the superclean limit it is given by \{9)=Gc{9)^ shown in Fig. 3. 

We can now calculate the thermal conductivity within a similar approximation for T<^ A 
and obtain 

«c 2 vl{eH) ^ 

for the superclean hmit (10) 



for the clean limit 

where Knc is the thermal conductivity in the normal state along the c axis while 
«;oc='^c(H=0) in the clean limit. The linear H dependence of the heat current is the signature 
of the superclean limit. In the clean limit, on the other hand, the thermal conductivity starts 
from the nonvanishing value for H=0. Also the H dependence is sublinear in H. Similarly 
the heat current parallel to the a axis is given by 

2vl{eH) 
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FIG. 3. Angular function l{9)=Gc{0)^ ~ g(^) which determines the specific heat in the 
superclean limit for a=l and 0.6 and different node line positions xo=ckz- 

for the superclean limit (11) 



for the clean limit 



where 



TT JO 

[cos^ 9 + sin^ 6'(sin^ + a sin^ xo) 

+y/asmxocos(j)sm{26)]2 (12) 
Fa{0) = (cos^ e + a sin^ 6*)^ [cos^ e + {^ + a sin^ xo) sin^ 0] 

The 9- dependence of Fa(^) and G„(^^)= I{6)I{0) for a=l and x=0, f and | are shown in 
Figs. 1 and 2 (full lines) and similar curves for a;=0.6 in Figs. 4 and 5. The boundary values 
for Xo=0 are Fa(|) = and Ga(|) = -^-Ja. From the superclean limit in both cases 
a and c in Figs. 2 and 5 we notice that the angular dependence is decidely nonmonotonic 
for the even gap function as opposed to the monotonic decrease of Gq c (Xo 7^0) in the odd 
parity case (xo=0). Therefore thermal conductivity measurements in UPd2Al3 should be 
able to decide definitely between the two cases. So far we have neglected the phononic 
thermal conductivity, which may be easily considered as in |Q. 

Comparison of a=l and a=0.6 results show that the difference in 9- dependent curves 
Fa,c(6') and Ga,c{0) for the various node positions xo becomes smaller with decreasing a. For 
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FIG. 4. Angular functions Fa{0) and Fc{0) in the clean limit and anisotropic case (a= 0.6) 
for different node line positions xo=c^z- 



specific compounds the value of a may be obtained from the upper critical field anisotropy 
according to Vc/va=Hc2||/Hc2±. For a small value a=0.1 the difference at 6' = | (for the 
various Xo) becomes also quite small (~10%). This means that in the very anisotropic 
cases like Sr2Ru04 {y/a = ^, i.e. a=0.0025) the angular 9 dependence therefore can only 
determine the proper (ab-) plane of node lines but not the position along the c- axis. On 
the other hand UPd2Al3 with a= 0.69 is a much more favorable case and the 6- dependence 
of c, Cs should render also the node position along the c- axis. 

In conclusion there are many experimental techniques to discover the nodal supercon- 
ductivity. In particular the a/^- dependence of the specific heat provides the clearest signal 
for it @J2^. On the other hand to determine the nodal directions in A(k) we need more deli- 
cate investigations like the phase sensitive experiments . However, the elegant tricrystal 
experiments appear to not transport to other nodal superconductors aside from high Tc 
superconductors. In this circumstance the angular dependent thermal conductivity will pro- 
vide a unique window to look at the nodal directions . From the angular dependence of 
the thermal conductivity Izawa et al have succeeded in deducing the symmetry of A(k) in 
Sr2Ru04 [Q, CeCoIns ^ and more recently k-(ET)2Cu(NCS)2 0- Therefore it is of great 
interest to explore this technique to still unidentified order parameter in UPd2Al3 ,UNi2Al3 
and other nodal superconductors. 
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FIG. 5. Angular functions Ga{0)= 1{9)I{9) and Gc{6)= 1(6*)^ in the superclean limit and 
anisotropic case (a=0.6) for different node line positions xo=ck^- 
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